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COMPUTER AIDED SOLUTION OF THE INVARIANCE EQUATION 
FOR TWO- VARIABLE STOLARSKY MEANS 

SZABOLCS BAJAK AND ZSOLT PALES 

Abstract. We solve the so-called invariance equation in the class of two-variable Sto- 
larsky means {S p . q : p,q € K.}, i.e., we find necessary and sufficient conditions on the 6 
parameters a, b, c, d,p, q such that the identity 

S P ,q(S a ,b{x,y) ) S Cjd (x,y)) = S p , q (x,y) (x,y e M+) 

be valid. We recall that, for pq(p — q) ^ and x ^ y, the Stolarsky mean S p q is defined 
by 

S pq{ x,y):=(^ P -y P \ 

In the proof first we approximate the Stolarsky mean and we use the computer algebra 
system Maple V Release 9 to compute the Taylor expansion of the approximation up to 
12th order, which enables us to describe all the cases of the equality. 



1. Introduction 

Let R + denote the set of positive real numbers throughout this paper. A two-variable 
continuous function M : R+ — > R + is called a mean on R + if 

(1) min(x, y) < M(x, y) < max(x, y) (x, 

holds. If both inequalities in ([1]) are strict whenever x ^ y, then M is called a strict mean 
on R + . 

Given three strict means M, N, K : R + , we say that the triple (M, N, K) satisfies 

the invariance equation if 

(2) K(M{x, y), N{x, y)) = K(x, y) (x, y e R+) 

holds. If 02]) holds then we say that K is invariant with respect to the mean-type mapping 
(M, N). It is well known that K is uniquely determined by M and N, and it is called the 
Gauss composition K = M <g> of M and N. For this terminology and result, see for 
instance the papers [15], |9]. 

The simplest example when the invariance equation holds is the well-known identity 



x + y 2xy 
V 2 x + y 

that is, 

S(.x, y) = 9(A(x, y), X(x, y)) (x, y G R+), 
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where A, 9, and "K stand for the two- variable arithmetic, geometric, and harmonic means, 
respectively. Another less trivial invariance equation is the identity 

A®S(x,y) — A® 3(A(x,y), 5(x, y)) (x,y e R+), 

where A ® 9 denotes Gauss' arithmetic- geometric mean defined by 



A®5(x,y) 



dt 



x 2 cos 2 1 + y 2 sin t 



[x,y G 



The reader is recommended to consult the book [3j for more details and history of this 
deep theory. 

The invariance equation in more general classes of means has recently been studied ex- 
tensively by many authors in various papers. The invariance of the arithmetic mean A (i.e., 
when in ([2]) K is the arithmetic mean) with respect to two quasi-arithmetic means was 
first investigated by Matkowski [15] under twice continuous differentiability assumptions 
concerning the generating functions of the quasi-arithmetic means. These regularity as- 
sumptions were weakened step-by-step by Daroczy, Maksa, and Pales in the papers [7], [8], 
and finally this problem was completely solved assuming only continuity of the unknown 
functions involved [9|. The invariance equation involving three weighted quasi-arithmetic 
means was studied by Burai [5J and Jarczyk-Matkowski [14] , Jarczyk [T3]. The final an- 
swer (where no additional regularity assumptions are required) has been obtained in |13| . 
In a recent paper, we have studied the invariance of the arithmetic mean with respect to 
two so-called generalized quasi-arithmetic mean under four times continuous differentia- 
bility assumptions |2]. The invariance of the arithmetic mean with respect to Lagrangian 
means was the subject of investigation of the paper |17| by Matkowski. The invariance of 
the arithmetic, geometric, and harmonic means with respect to the so-called Beckenbach- 
Gini means was studied by Matkowski in [16] . Pairs of Stolarsky means for which the 
geometric mean is invariant were determined by Blasiriska-Lesk-Glazowska-Matkowski 
[6]. The invariance of the arithmetic mean with respect to further means was studied by 
Glazowska-Jarczyk-Matkowski |12| . Burai jl] and Domsta-Matkowski [10J. 

An important class of two variable homogeneous means are the so-called Stolarsky means 
(cf. Stolarsky [IS]). Given two parameters p, q G M, the two- variable mean S Piq : IR+ — > M+ 
is defined by the following formula 



S p ,q{x, y) :- 



q(x p — y p ) 
p(x q — y q ) 



i 

P-9 



exp 



1 

P 



x p log x — y p log y 
x p — y p 



p(\ogx 
g(logx 



logy) 
log?/) 



if (p - q)pq ^0,x^t/ 

if V = 1 7^ 0, x ^ y, 
if V ^ 0, q = 0, x ^ y, 



x q — y q 



if p ■ 

if p ■ 
if x 



0, q ^ 0, x ^ y, 
q = 0, 

y, 



for every x,y E E + . 

It can easily be seen that the power (or Holder) mean of exponent p can be obtained 
as S2 P ,p- In particular, Sq^ and S'_2._i are the arithmetic, geometric, and harmonic 
means, respectively. Moreover, it follows by a simple computation that S^q is equal to the 
geometric mean whenever p + q = 0. One can also check that the Stolarsky means are 
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symmetric with respect to their variables and their parameters as well, i.e., S P) q(x,y) = 
S Pi q(y, x) = S q>p (x, y) for all x, y G R + and p, q G R. 

The aim of this paper is to solve the invariance equation in the class of Stolarsky means, 
i.e., to solve ([2]) when each of the means M, N and K is a Stolarsky mean. More precisely, 
we want to describe the set of all 6-tuples (a, b, c, d,p, q) such that the identity 

(3) Sp !q (S a>b (x, y), S c4 (x, y)) = S P)i (x, y) (x, y G R+) 

holds. The analogous problem concerning the so-called Gini means, i.e., the description of 
the set of all 6-tuples (a, b, c, d, p, q) such that 

(4) G Pjq (G a>b (x,y),G Cid (x,y)) = G m (x,y) (x,y eR+) 

be valid, has recently been solved in the paper [1] by the authors. Recall that given two 
parameters p, q G R, the two- variable mean G p>q : R+ — > R + is defined by the following 
formula (cf. Gini 



G P Jx,y) 



X P + yP 
x q +y q 



1 

p-q 



for p 7^ q, 



( x p log x + y p log y \ 

ex P V7~v for p = q, 

\ xp + yP J 

for x, y G R+. The class of Gini means is also a generalization of the class of power means, 
since taking q = 0, we immediately get the power (or Holder) mean of exponent p. The 
main result of the paper [lj is contained in the following theorem: 

Theorem G. Let a,b,c,d,p,q G R. Then the invariance equation (J3J) is satisfied if and 
only if one of the following possibilities hold: 

(i) a + b = c + d = p + q = 0, i.e., all the three means are equal to the geometric mean, 

(ii) {a, b} = {c, d} = {p, q}, i.e., all the three means are equal to each other, 

(Hi) {a, b} = {—c,—d} andp + q = 0, i.e., G ps is the geometric mean andG a ^ = G_ Cj _d ; 
(iv) there exist u, v G R such that {a, b} = {u + v, v}, {c, d} = {u — v, —v}, and {p, q} = 

{u,0} (in this case, G PtQ is a power mean), 
(v) there exists w G R such that {a, b} = {3w,w}, c + d = 0, and {p, q} = {2w,0} (in 

this case, G p>q is a power mean and G Cy d is the geometric mean), 
(vi) there exists w G R such that a + b = 0, {c,d} = {3w,w}, and {p,q} = {2w,0} (in 

this case, G p>q is a power mean and G a ,b is the geometric mean). 

As an obvious consequence, we obtain the following solution for the so-called Matkowski- 
Suto equation, i.e., when G p>q is equal to the arithmetic mean in the invariance equation 
01]), which happens exactly when {p, q} = {1,0}. 

Corollary G. Let a, b, c, d G R. Then the Matkowski-Suto equation 

G atb (x,y) + G Ctd (x,y) = x + y (x,yGR + ) 

is satisfied if and only if one of the following possibilities hold: 

(i) {a, b} = {c,d} = {1,0}, i.e., the two means are equal to the arithmetic mean, 
(ii) there exists v G R such that {a, b} = {1 + v, v}, {c, d} = {1 — v, —v}, 
(Hi) {a, b} = {|, |} and c + d = (in this case, G c ^ is the geometric mean), 
(iv) a + b = and {c, d} = {|, |} (in this case, G a ^ is the geometric mean). 

The approach followed in [lj heavily used the Maple computer-algebra system to perform 
the tedious computations of the various partial derivatives of Gini means up to 12th order. 

The main result of this paper completely solves the invariance equation in the class of 
Stolarsky means. 
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Theorem S. Let a,b,c,d,p,q G R. Then the invariance equation ([3]) is valid if and only 
if one of the following possibilities holds: 

(i) a + b = c + d = p + q = 0, i.e., all the three means are equal to the geometric mean, 
(ii) {a, b} = {c,d} = {p, q], i.e., all the three means are equal to each other, 
(Hi) {a, b} = {— c, — d} and p + q = 0, i.e., S Pjq is the geometric mean and S a ^ = S- c -d- 

As a consequence, we obtain the following solution for the Matkowski-Suto equation, 
i.e., when S p>q is equal to the arithmetic mean in the invariance equation ([3]), which happens 
exactly when {p, q] = {2, 1}. 

Corollary S. Let a, b,c,d G R. Then the Matkowski-Suto equation 

S a>b (x,y)+S Cjd (x,y) = x + y (x,yGR+) 

holds if and only if {a, b} = {c,d} = {2, 1}, i.e., both means are equal to the arithmetic 
mean. 

It is interesting to observe here that the parameter sets when (jl]) holds is much bigger 
then the corresponding set for Q). 

In the proof of the above results first we construct an approximation of the Stolarsky 
means. Then we use the computer algebra system Maple V Release 9 to compute the 
Taylor expansion of a function (which is in terms of the approximated means) up to 12th 
order. Finally, this enables us to describe all the cases of the equality. 



2. The proof of Theorem S 
First we recall the characterization of the equality of two variable Stolarsky means. 

Lemma 1. (Cf. [IS] ) Let a, fe, c, d G R. Then the identity 

(5) S a ,b(x, y) = S c4 (x, y) (x, y G R+) 

holds if and only if one of the following possibilities is valid: 

(i) a + b = c + d = and, in this case, the two means are equal to the geometric mean, 
(ii) {a, b} = {c, d}. 

The calculations of the higher-order partial derivatives of Stolarsky means at the point 
(x,y) = (1,1) is too complicated (even for computer-algebra systems) because the main 
expression that defines these means is singular at diagonal points. Therefore, to simplify 
these computations, we will approximate the Stolarsky mean up to a sufficiently high order. 
First, we express the Stolarsky mean in another form. 

Lemma 2. Forp,q £ R and x,y G R + , the Stolarsky mean S pq (x,y) can be rewritten as 
(6) 

s 1 

L(p log x,p logy) 



S p ,q(x, y) 



L(qlogx,q\ogy) 



ifP 9, 



. diL(p log x,p logy) log x + d 2 L(p hgx, p logy) logy \ 

ex P TT~\ \ \ VP = 9' 

Lyp log x,p logy) 



where the function L : R — >■ R is defined by 



CO 



(7) L(u,v):=J2 



u n-l + u n~2 v ^ 1_ uv n-2 + 



n ! 

n=l 
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Proof. Utilizing the Taylor series expansion and substituting u := logx and v := logy, we 
have, for p ^ and x ^ y, 



X P _ yP 



p (log x — log y) p(u — v) p{u — v) 

1 



pu (pu) 2 _ 1 _ pv _ (pv) 2 
1! 2! 1! 2! 



u — V 



u — v p(u 2 — v 2 ) p 2 (u 3 — v 3 ) 

+ ^-r, + 7T. ~ + 



1 ! 



2! 



3! 



p(u + v) p 2 (u 2 + uv + v 2 ) 

= 1 + — + — + ■ ■ 

2! 3! 

Similarly, we can get that L(pu,pu) = e pu , hence 



X P _ yP 



L(pu,pv). 



if p ^ 0,x ^ y, 



L(plogx,plogy) := < 



x p , 



p (log x — log y) ' 

if a; = y, 
if p = 0. 

Using the above identity and the definition of the Stolarsky means, equality ([6]) follows 
immediately for p ^ q. The formula for the case p = q is derived by taking the limit 
q — > p. □ 

In order to obtain high-order approximation of Stolarsky means, for fceN, define 



(8) 



L k (u,v 



n=l 



n \ 



(u, v G 



and, for x, y G R + , p, g G 
(9) 

✓ 

'L fc (plogx,plog?/)^ 

Sj, g (a:,3/) := < 



cxp 



L k (q\ogx,q\ogy) 

d 1 L k (p\ogx,p\ogy)\ogx + <9 2 L fe (p log x,p logy) logy 



if p ^ g, 
if p = g. 



L fc (p log x,p logy) 
Lemma 3. For p, q G R, G N, 

(10) 9^^(1,1) = 9^^(1,1) 

« > 1 , j > 1 and z + j < fc. 

Proof. Observe that, for a > 0, > with a + j3 < k, we have 

d?fl£L(0,0) = d?dlL k (0,0). 

In the case p 7^ q, by applying the elementary rules of differentiation and the identities 
and (P), we get that the formulae for the partial derivatives d\d^ 2 S Ptq (l, 1) and d^d^SpJl, 1) 

are the same expressions of all the partial derivatives d^d^ 1/(0,0) and di^Lk(0,0), re- 
spectively, where < a < i and < < j. Hence, equality (TTOl) follows for p 7^ g. In the 
case p = g, the identity is derived by taking the limit q — » p. □ 

To make the calculations even more simple, we define the function E : R 4 — > R, by 

' log L(pu, pv) — log L(qu, qv) 



E(p,q,u,v) :=log {S p , q (e» , e v )) 



p-q 

d\L(pu,pv)u + d2L(pu,pv)v 
L(pu,pv) 

The invariance equation ([3]) can be rewritten in the form 

E(p, q, E(a, b, u, v), E(c, d, u, v)) = E(p, q, u, v) (u, v G 



if P ^ q, 
if p = q. 
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hence 

F(x) := E(p, q, E(a, b, x, —x), E(c, d, x, —x)) — E(p, q, x, —x) = (x 6 1). 

Thus, the derivatives F^ vanish at x = for all m > 0, m £ 7L. By the symmetry of 
Stolarsky means, F is an even function. Therefore, F^ 2m+1 \0) = holds automatically. In 
order to obtain the necessity of the conditions of the theorem, we will need to investigate 
the equalities 

(11) F (2m) (0) = 0, (l<m<6). 

To compute these derivatives, define the function F k : K — > R, k G N, as 

F k (x) := E k (p,q,E k (a,b,x,-x),E k (c,d,x,-x)) - E k (p, q, x, -x), 

where 

' log L k (pu,pv) - log L k (qu,qv) 



E k (p,q,u,v) :=log {S k m (e u ,e v )) 



p-q 

diL k (pu,pv)u + d 2 L k (pu,pv)v 



if p ^ q, 
if p = q. 



L k (pu,pv) 

It follows from Lemma [3] that the derivatives of the functions E k and E at (p, q, 0, 0), and 
hence the derivatives of F k and F up to the order k at the point x = coincide. Therefore, 
to analyze the equations (ITT]) , it is sufficient to consider the identities 

(12) if m) (0) = 0, (l<m<6), 

where 2m < k. Thus in the proof of the theorem, we compute the Taylor expansion of 
the appropriate F k at x — 0, hence we can obtain conditions for the unknown parameters 
a, b, c, d,p, q. 

The proof of Theorem S. Assume that ([3]) holds. In the syntax of the Maple language, we 
define the functions L k , E k and F k in the following way: 



=(u,v J k)->add((l/n!)*add(u-(n-i)*v-(i-l) ,i=l. .n) ,n=l. .k) ; 
=(p,q,u,v,k)->(l/(p-q))*(ln(L(pu,pv,k))-ln(L(qu,qv J k))) ; 
=(x,k)->E(p,q,E(a,b,x,-x,k) ,E(c ,d,x, -x,k) ,k) -E(p, q,x, -x,k) ; 



This produces the following output: 

L :=(u, v, k) — > add ^ p- -, n = 1 . . k 

log (L(pu,pv } k)) - log (L(qu, qv, k)) 
p-q 

F :—(x, k) — > E(p, q, E(a, b, x, —x, k), E(c, d, x, —x, k),k) — E(p, q, x, —x, k). 
First we evaluate the second-order Taylor coefficient C2 of F at x = by replacing F 

i».v /•:•',: 

> C[2] :=simplify(coeftayl(F(x,3) ,x=0,2)) ; 

1 j 1 1 1 1 1 
Go = — d q -\ p q H a H c. 

12 6 q Yl y 6 y 12 12 
(Note that the Maple-definition of the function E is valid only if (p — q)(a — b)(c — d) ^ 0, 
however, the Taylor coefficient C2 and also the subsequent ones, are correct also in the 
singular case (p — q)(a — b)(c — d) = 0.) 

By the invariance equation C2 = 0, therefore 

a + b + c + d p + q 
4 ~ 2 



INVARIANCE EQUATION FOR TWO- VARIABLE STOLARSKY MEANS 7 

In order to simplify the evaluation of the higher-order Taylor coefficients, we introduce the 
notations 

a + b + c + d p + q 



w 



v 



4 

a + b — (c + d) 



4 

/ n — r 

(13) t : 



p — q\ 2 



2 

(a - bf + (c - df 

r : = , 

8 

(a-bf -(c-d) 2 

S '= 8 • 

(In the definition of w we utilized the condition C2 = 0.) Then we can express the 
parameters a, b, c, d, p, q in the following form: 

> a: =w+v+sqrt (r+s) ; b : =w+v-sqrt (r+s) ; 
c : =w-v+sqrt (r-s) ; d : =w-v-sqrt (r-s) ; 
p : =w+sqrt (t) ; q: =w-sqrt (t) ; 



ir t v + y/r + s 



w + v — y/r + 



w — v + y/r 



w — v — y/r — s 
p := w + y/t 
q := w — yft 

Now we evaluate the 4th order Taylor coefficient C4 of F (i.e., of F5) at x = by inputting: 
> C[4] :=simplify(coeftayl(F(x,5) ,x=0,4)) ; 

1 4 2 1 1 

C4 := -—vs - ——wv + — wt - —wr 
45 135 45 45 

4 

The condition C4 = yields that wt = wr + v s H — wv 2 . 

3 

If w — 0, then p + q = 0, which means that S P) g is the geometric mean. Therefore, the 
invariance equation can be written as 

Sa,b(x, y)S c , d (x, y) = xy (x, y E K+). 

This results 

S a ,b{x, y) = . ) . r = S- c ,- d (x, y) (x, y e R+). 

Sc,d{i/x, yy) 

Using Lemma [U this identity yields that either a + b = c + d = 0oi {a, b} = {— c, —d} 
must hold. In this case we get that one of the conditions (i) or (iii) of our theorem is valid. 
Conversely, if conditions (i) or (iii) hold then (J3]) can easily be seen. 

In the rest of the proof, we may assume that w is not zero. Then, from condition C4 = 0, 
we can express t in terms of w,v,r,s: 



> t :=r+v*s/w+(4/3)*v~2; 

(14) t:=r + — + - 

w 3 

Next, we evaluate the 6th order Taylor coefficient C$ of F (i.e., of F 7 ) at x = 0: 



.2 



v 
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> C[6] :=simplify(coeftayl(F(x,7) ,x=0,6)) ; 



We get that 

_ 2(9w 2 s 2 + 8w 2 v 4 + 45w 2 rv 2 + 39w 3 sv + 6wv 3 s - 13w 4 v 2 - 9v 2 s 2 ) 
6 "~ 8505w 

If v = 0, then Cq = implies that s = 0. Hence, from ( Tl4|) it follows that t = r and we 
get that a = c = p and b = d = q, i.e., condition (ii) of our theorem holds. Conversely, if 
condition (ii) holds then the invariance equation ([3} is trivially valid. 

In the rest of the proof, we may assume that v is also not zero. Observe that the 6th 
order coefficient Cq does not involve higher-order powers of r. Therefore, the equation 
Cq = can be solved for r. Temporarily, we denote this solution by R: 



> R: = (13*w"4*v'-2-9*w~2*s"2-8*w~2v"4+9*v'~2*s~2-39*w'~3*v*s-6*w*v"3*s)/ 
(45*w~2*v~2) ; 



„ 13w 4 v 2 — 9w 2 s 2 — 8w 2 v 4 + 9v 2 s 2 — 39w 3 vs — 6wv 3 s 
( 15 ) R: = az 9 9 

Finally, we evaluate the 13th order Taylor polynomial of F 13 at x = (the Maple output 
is suppressed by putting : instead of ; to the end of the Maple command, for the sake of 
brevity), then we extract the 8th, 10th and 12th order Taylor coefficients, denoted by Cg, 
Cio and C12, respectively, and replace r by R by inputting: 

> T:=simplify(taylor(F(x, 13) ,x=0,13)) : 
for i from 8 to 12 by 2 

do C [i] : =simplif y(subs (r=R, simplif y(coef f (T,x, i) )) ,f actor) od; 



C 8 := — — (4347wV - 3616wV + 1242wVs 2 + 16137iuW + 224wV 

95681 25ur%^ v 

- 4QUw 3 v 3 s 3 + 75Qw 3 v 7 s + 81wV + 11976™ Vs + 4992wV - lQ2w 2 v 2 s 4 

- 17379wVs 2 - 12732wVs + 297wv 5 s 3 + 8lv 4 s 4 ) 



C w := — (141632m; 1 V - 2187wVs 6 - 58806wVs 5 + 107406wVs 5 

2841733125w; 5 f 4 v 

- 338432w 1 V + +34928wV 2 - 523908wVs 3 + 2187wVs 6 + 126459wW 
+ 316272tuV° + 143541wVs 2 + 92508wV 1 s + 24948wVs 4 + 131901Gw n v B s 

- 847044wVs + 927324w 1 %V - 1197324wVs 2 - 564480wVs - 729w 6 s 6 
+ 55026w Vs 3 - 794610^ Vs 3 + 672786wVs 4 + 3402ttw V - 52002w W 

- 313065wVs 4 + 1263492w;Vs 3 - 384669wVs 4 + 729w 6 s 6 ) 
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C l2 := — (54403812wVs 8 - 36269208wVs 8 - 2472351012m W 

872767286015625wV V 

- 36269208wVs 8 - 6728400999u>Vs 6 - 23518313469u?Vs 6 + 28808296272u> 1 V 4 

- 43502512128w 1 V - 32563674368w 1 V° + 55532868672u4V 2 + 408258 1552u?V 6 
+ 51381378WV + 1616450580w W + 66139190820wVs 5 + 5633126127wW 

- 27686800350u>Vs 5 + 23171026860W 1 Vs 5 - 198268186320u> 1 Vs 3 + 9067302w 8 s 8 

- 50400222003u> 1 Vs 4 + 8922828051u> 1( W + 527871816™ W + 2575113768iuVi 
+ 33853773582w 1 Vs 4 + 8610065640w W + 78734957472wVs 4 + 806989878w 9 t;s 

- 56327456871™ W - 297179050266w; 1 V s 2 - 14684723544w 1 Vs 
+ 203758408416w 1 Vs - 195782704164u4 V x s - 67821635178™ W 
+ 10150547496wV 5 s + 13660568109™ W - 84844566450™V 1 s 3 
+ 337457994084™ 1 Vs 2 + 2387944944w 1 Vs 2 - 3441528204™ V 3 s 

- 67735201590™Vs 3 + 342237888720™ 1 Vs 3 - 961134012™Vs 7 
+ 2079601460l™Vs 6 - 63239867910™Vs 5 + 9067302v 8 s 8 ) 

The Taylor coefficients C%, C w , and C± 2 are of the form 

12 2 
n — p n — p n _ p 

8 9568125™ 3 *; 2 8 ' 10 2841733125™ 5 *; 4 10 ' 12 8727672860 15625™ 7 *; 6 12 ' 

where P 8 , P 10 , and P 12 are polynomials of the variables v,w, s. They can be obtained by 
the following Maple commands (whose output is suppressed): 

> P[8] :=op(2,C[8]) : P [10] :=op(2,C [10] ) : P [12] :=op(2,C [12] ) : 

The equalities C 8 = C w = C\ 2 = imply that P 8 = Pio = P\ 2 = 0. 

The variable s is a common root of the polynomials P 8 and P\q. Therefore the resultant 
#8,10 °f these two polynomials (with respect to s) is zero: 

> R [8 , 10] : =f actor (resultant (P [8] , P [10] , s) ) ; 

i? 8 ,io := 28242953648100000000w 2 V 4 (t; - wf{v + wf (395726752304 v 32 - 28019198519832 w 2 v 3 
+ 1192972799035666 w 4 v 28 - 36617671790074251 w 6 v 26 + 601554420387156651 w 8 v 24 

- 3652037976710860175 w 10 v 22 - 1101310194408221307 w 12 v 20 + 62048533824813847173 w 14 v 18 

- 175575191501013599783 w 16 v 16 + 52614376847529172973 w 18 v 14 + 435211540238087039223 w 20 

- 793895884964266327270 w 22 v w - 773252618095825970136 w 24 v 8 - 492199682627262911866 w 26 
+ 183522699320559043726 w 28 v 4 - 43030934088053846752 w 30 v 2 + 912066926976343384 w 32 ) 

The resultant is zero if either vw(v — w)(v + w) — holds or v and w are solutions of a 
homogeneous two variable polynomial equation of degree 32. Writing w in the form 

> w:=z*v; 



w := zv 

we get that 2 is a root of a 32nd degree polynomial #3,10, where: 
> P[8,10] :=simplify(op(4,R[8,10])/v~32); 
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P 8 ,io := 395726752304- 28019198519832/ + 1192972799035666/ - 36617671790074251/ 
+ 601554420387156651/ - 3652037976710860175/ - 1101310194408221307/ 2 
+ 62048533824813847173/ 4 - 175575191501013599783/ 6 + 52614376847529172973/ 8 

- 435211540238087039223/ - 793895884964266327270/ 2 + 773252618095825970136/ 4 

- 492199682627262911866/ 6 + 183522699320559043726/ 8 - 43030934088053846752/ 

+ 912066926976343384/ 2 . 

The variable s is also a common root of the two polynomials Pg and Pu- Therefore the 
resultant i?8,i2 of these polynomials (with respect to s) is again zero: 

> R[8,12] :=f actor (resultant (P [8] ,P[12] ,s)) : 

We now get that v and w are solutions of a homogeneous two variable polynomial of degree 
44, whence we get that z is a root of the 44th degree polynomial Ps,i2, where: 

> P[8,12] :=op(4,R[8,12]) ; 

Now computing the resultant of the two polynomials Ps,io an d Ps, 12 by 

> Q : =resultant (P [8 , 10] , P [8 , 12] , z) : 

it follows that Q is a (huge) nonzero number, hence Ps,io an d Ps,i2 cannot have a common 

root. This proves that Rs,xo and -Rs,i2 can be simultaneously zero if and only if vw(v — 

w)(v + w) = holds. Hence Cg = C\q = C\i = can hold only in this case. 

Since we have vw 7^ 0, hence (v — w) (v + w) — must hold, i.e, v = ±w. Thus, from 

w 2 13 
(Il5p and f )14p . we get that r = — =)= s and t = — w 2 , respectively. 

9 9 

In the case when v = w, the equations in (|T3l) yield that 

(16) c=-d, a = ^w, b= ^ w > P= ( 1 + _ 3~) W ' g= ( 1_ "T^) U '' 

The first equality yields that S c ,d is the geometric mean. Hence, we may assume that 
c = —d = 1. To simplify the computations, we can also assume that w — 3. We show that 
these parameters are not solutions of the invariance equation. For k G N, we now have 
that 

F k (x) := E k (3+V13, 3- Vl3, E k {7, 5, x, -x), E k {l, -l,x, -x)) -E k 3-Vl3, x, -x) . 

In Maple, we input 

> F:=(x,k)->E(3+sqrt(13) ,3-sqrt(13) ,E(7, 5,x, -x,k) ,E(1 , -1 ,x, -x,k) ,k) - 
E(3+sqrt(13) ,3-sqrt(13) ,x,-x,k) ; 

We compute the 10th order Taylor coefficient of Fn by 

> simplify (coeftayl(F(x, 11) ,x=0,10)) ; 

whence we get that this coefficient is — — , i.e., it is not zero, which means that the 

5775 

parameters in (TIB"]) do not provide solution to the invariance equation. 
In the case when v = —w, we have that 



7 ,5 / V13 / 1 V13\ 
a = -b, C= 3 W ' 3 W ' p= \ ^r) W ' <?=|1 g- I w, 

whence a similar calculation as in the previous case shows that we again do not get an 
additional solution to the invariance equation. □ 
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